Could the next generation of cosmology experiments exclude supergravity ? 
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Gravitinos are expected to be produced in any local supersymmetric model. Using their 
abundance prediction as a function of the reheating energy scale, it is argued that the next 
generation of Cosmic Microwave Background experiments could exclude supergravity or strongly 
favor "thermal-like" inflation models if B mode polarized radiation were detected. Galactic 
cosmic-ray production by evaporating primordial black holes is also investigated as a way of 
constraining the Hubble mass at the end of inflation. Subsequent limits on the gravitino mass and 
on the related grand unification parameters are derived. 
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INTRODUCTION: GRAVITINOS IN THE 
EARLY UNIVERSE 



Although not yet experimentally discovered, super- 
symmetry (SUSY) is still the best - if not the only 
- natural extension of the standard model of particle 
physics. It could provide a general framework to 
understand the origin of the fundamental difference 
between fermions and bosons and could help to resolve 
the difficult problem of mass hierarchies, namely the 
instability of the electroweak scale with respect to 
radiative corrections. In global supersymmetry, the 
generator spinors f are assumed to obey = 

If one wants to deal with local supersymmetry, 
or supergravity, this condition must be relaxed and ^ 
becomes a function of the space coordinates x. New 
terms, proportional to d^^{x), must be canceled by 
introducing a spin 3/2 particle, called gravitino, as 
vector bosons are introduced in gauge theories. The 
gravitino is part of an N=l multiplet which contains 
the spin 2 graviton (see for an introductive review) 
and, in the broken phase of supergravity, super-Higgs 
effects make it massive through the absorption of the 
Nambu-Goldstone fermion associated with the SUSY 
breaking sector. 

It has long been known that if the gravitino is 
unstable some severe constraints on its mass must be 
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considered in order to avoid entropy overproduction j^]: 
1TT-3/2 <L 10 TeV. On the other hand, if the gravitino is 
stable, its mass should satisfy < 1 keV to keep 

the gravitinos density smaller than the full Universe 
density (^13/2 < ^tot)- In spite of the huge dilution, 
those constraints are not fully evaded by inflation as 
gravitinos should be reproduced by scattering processes 
off the thermal radiation after the Universe has reheated 
H H S H H m |l3- As the number density of such 
secondary gravitinos is expected to be proportional to 
the reheating temperature, it is possible to relate the 
energy scale of inflation with the requirement that they 
are not overproduced. 



In the ffist part of this paper, the next generation of 
cosmic microwave background (CMB) detection experi- 
ments is considered as a way of possibly excluding su- 
pergravity. It is shown that the energy scale of inflation 
required to produce an observable amount of tensor mode 
in the background radiation is not compatible with local 
supersymmetry in the standard cosmological scenario. In 
the second part, a new way of constraining the grav- 
itino mass, based on evaporating primordial black holes, 
is investigated. Taking into account that the black hole 
masses cannot be much smaller than the Hubble mass at 
the formation epoch, it is suggested that a detection of 
cosmic-rays produced by the Hawking mechanism would 
lead to a lower bound on the reheating scale and, there- 
fore, on the gravitino mass. Links with grand-unified 
models are given, as an example, in the conclusion. Fi- 
nally, the basics of the propagation model used to relate 
the source term to the local spectrum are given in the 
Appendix A. 
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II. TENSOR MODE IN THE COSMOLOGICAL 
BACKGROUND 

Observational cosmology has recently entered a new 
era thanks to several experiments dedicated to the CMB 
measurements (see [H^), e.g. Maxima, BOOMERanG, 
ACBAR, DASI, CBI, VSA, ARCHEOPS, and WMAP. 
They give strong evidences in favor of the inflationary 
scenario: a density extremely close to the critical value, 
a nearly scale invariant power spectrum, and a gaussian 
structure of the perturbations. Furthermore, in addi- 
tion to the temperature anisotropics, the polarization of 
the CMB has also been recently observed 0, . For 
the time being, only the even-parity E mode has been 
detected and the odd-parity B mode is still to be dis- 
covered. The latter is of specific importance as it would 
probe the primordial gravitational waves through tensor 
perturbations. Their amplitude can be expressed with 
the Hubble parameter and the potential of the scalar field 
driving inflation |T^ : 
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The amplitude of the polarization B mode is therefore 
directly proportional to Ein/i- 

Figure ^ shows the 1 a sensitivity of the Planck 
satellite to polarization, as computed with CMBfast 
On the same plot, the B mode polarization in a 
standard ACDM cosmology with an inflation energy 
scale Einfi ~ lO^^GeV (dotted line) is also represented. 
Increasing (resp. lowering) Einfi would result in increas- 
ing (lowering) the amplitude of the primordial B mode 
thus making it easier (more difficult if not impossible) 
to detect. On the contrary, the level of the expected 
B mode induced by weak lensing is fixed and rather 
accurately predicted since it results from lensing effects 
on the polarization E mode due to scalar perturbations 
that are now well constrained. The challenge in the 
detection of the primordial B mode and the estimation 
of Einfi is then to have a sensitive enough experiment 
and to avoid contamination by weak lensing. For the 
Planck experiment, the major hope is the detection 
at low i thanks to the high reionization optical depth 
suggested by WMAP 16] . In the case of limited sky 
coverage experiments, the weak lensing contribution will 
have to be removed. 

With the Planck sensitivity, the B mode should be 
detected (3ct) if E.nfi > 10^^ GeV [H [13. This case 
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FIG. 1: Sensitivity (1 a) to polarization of the Planck satel- 
lite (solid line) versus the expected B mode polarization in a 
standard ACDM cosmology with an inflation energy scale of 
^ W^^ GeV (dotted line). Planck should provide significant 
detection of this tensor mode, especially at low multipole £ 
where reionization boosts the power spectrum. The B mode 
induced by weak lensing is also represented (dot-dash line) 
and dominated the primordial spectrum for £ < 200. 



would be in severe conflict with most supersymmetric 
models. Indeed, in mSUGRA, the gravitino mass is, 
by construction, expected to lie around the electroweak 
scale, i.e., in the 100 GeV - 1 TeV range Consider- 
ing that Deuterium and "^He should not be overproduced 
by photodissociation of ^He below 700 GeV and that 
Deuterium should not be destroyed beyond the al- 
lowed observational values [23| above 700 GeV [2ll |. 
the reheating temperature must remain lower than 
2 X 10^ GeV if the branching ratio of gravitinos into 
photons and photinos is assumed to be unity and lower 
than 5 x 10^^ GeV with a conservative branching ratio of 
1/10. The large difference between those limits and the 
energy scale required to produce a measurable amount 
B mode polarization makes the exact value of the 
branching ratio of gravitinos into photons and photinos 
irrelevant. A detection of the polarization B mode by 
the Planck satellite would therefore disfavor mSUGRA 
in standard cosmology. 

In gauge-mediated SUSY breaking alternative sce- 
narios, mostly interesting in accounting for a natural 
suppression of the rate of flavor-changing neutral-current 
due to the low energy scale, the situation is even more 
constrained. In this case, gravitinos are the lightest 
supersymmetric particles and requiring their density not 
exceed the total density imposes an upper limit on Trh 
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between 10^ and 10^ GeV for masses between 10 MeV 
and 100 keV Alth oug h some refined models can 

relax those constraints [23 j local supergravity would, 
in this case also, be in serious trouble if the reheating 
temperature was high enough to be probed by the 
Planck experiment. 

A possible way to get around these conclusions is 
to assume that a substantial amount of entropy was 
released after the gravitinos and moduli production, that 
would dilute them according to the entropy conservation 
(n/s ~ cte). Such a scenario can be realized while 
keeping the inflationary scale high, e.g.^ in thermal 
inflation [3 H^- Some studies jla] even show that a 
wide modulus mass region (m$ « 10 eV - lO'' GeV) 
would be allowed but it requires in most cases a very 
small reheating temperature. Recently, the curvaton 
scenario has also attracted a considerable interest 
as it generates a huge amount of entropy through a 
scalar field that dominates the radiation at a given 
epoch. One can then argue that a detection of tensor 
mode polarization, would strongly favor "thermal like" 
inflation scenarios if supergravity is to remain as the 
preferred extension of the standard model of particle 
physics. Interestingly, if evidence in favor of local 
supsersymmetry were obtained either by colliders or by 
independent astroparticle experiments, this could even 
be a very promising observational signature for thermal 
inflation. 




FIG. 2: Possible reheating temperatures Trh as a function 
of the PBH density (normalized to the critical density) for 
different antideuteron flux at 100 MeV : 2 x 10"'^, 2 x 10"*, 
2 X 10-^ 2 X 10"^" from right to left in m-^s-^sr-^GeV"^ 



Fortunately, the Planck satellite is not expected to be 
the ultimate experiment to study the CMB polarization 
and several improvements can be expected in the future. 
However, as pointed out in Refs. [28L l29l |. there remains 
a lower limit to the removal of the polarization B mode 
foreground induced by gravitational Icnsing which sets at 
present time the lower limit on the detectable inflation 
scale to a few times 10"'^^ GeV. This scale remains, how- 
ever, particularly interesting if the fundamental scalars 
driving the phenomenon are related with grand unifi- 
cation since it lies around the GUT energy (between 
10^^ GeV and 3 x 10^^ GeV depending on whether su- 
persymmetry is considered or not). It therefore makes 
sense to improve the polarization sensitivity to reach the 
capability to probe the typical GUT scale where infiation 
could have occurred if the gravitino limit is ignored. 



III. COSMIC-RAYS FROM EVAPORATING 
BLACK HOLES 

Another interesting way to experimentally probe the 
reheating temperature would be to look for evaporating 
primordial black holes (PBHs). Such black holes should 
have formed in the early Universe if the density contrast 
was high enough on small scales. Many different possible 
scenarios have been suggested to allow for an important 
PBH density (see for a review): a dust-like stage 
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FIG. 3: Lower limit on the reheating temperature Trh as a 
function of the 100 MeV antideuteron flux. 



4 



[3l|, general first order phase transitions |33|, a scale 
in the power spectrum |33l l34l |. to mention only the 
currently most discussed possibilities. Such PBHs 
of mass M should evaporate, following a Plank-like 
spectrum with temperature T = hc^ / {IGnkGAI) , which 
was derived by Hawking |85l | using the usual quantum 
mechanical wave equation for a collapsing object with a 
postcollapse classical curved metric instead of a precol- 
lapse Minkowsky one. If those black holes are present in 
our galaxy (even with densities as low as ^pbh ~ 10~^), 
the emitted quanta should contribute to the observed 
cosmic-rays. Among them, two kinds of particles are 
especially interesting: antiprotons and gamma-rays. 
Antiprotons are useful because the astrophysical back- 
ground coming from spallation of cosmic-rays on the 
interstellar medium (so-called secondary particles) is 
very small (the ratio p/p is smaller than lO"'^ whatever 
the considered energy) and very well known 's^ . A tiny 
excess due to evaporating black holes could therefore be 
easily probed in the low energy range ^3 since the shape 
of the PBH spectrum is dominated by fragmentation 
processes and is then softer than the secondary spec- 
trum. Gamma-rays, coming both from direct emission 
and from the decay of neutral pions, take advantage of 
the very small optical depth of the Universe for ~ 100 
MeV radiation 'SS"]: the source emission can be probed 
up to redshifts z ^ 700. Furthermore, the signal to noise 
ratio is optimal at this energy as the PBH spectrum 
becomes softer {dN/dE oc E'^ dN/dE ex E'^) 
above 100 MeV (roughly corresponding to the QCD 
confinement scale) because of partons hadronization and 
integrated redshift effects [s^ . 

Using those cosmic-rays, the experiments are cur- 
rently sensitive to PBHs with masses between 10^^ and 
10^^ g. Those values can be intuitively understood as 
resulting from two opposite effects. On the one hand, 
the temperature favors the light (i.e. hot) black holes 
but their number density is very small : by integrating 
the Hawking flux over energy, it is straightforward to 
show that the mass spectrum must be proportional to 
below = 5 X 10^** g (the initial mass of a black 
hole whose lifetime is equal to the age of the Universe) 
whatever the details of the formation mechanism [40l| . 
This is mostly due to the fact that the low-mass be- 
havior is fully governed by the evaporation process, as 
obtained by writing dn/dM = (dn/dMi) x (dMJdM) 
where M stands for the current mass value and Mi 
for the initial one. The evolution term dMi/dM is 
simply determined from Mi « (3at -I- M^Y^^ where 
a « {7.8(is=i/2 + 3.1ds=i} x 10^** g^s~^ accounts for the 
number of available degrees of freedom with ds^i/2 = 90 
and ds=i — 27 in the standard model the other 

hand, the "number density" effect favors the heavy black 
holes but their low temperature makes the emission 
rate very small, especially when heavy hadrons are 
considered. 



The important point for this study is that only black 
holes formed after inflation would contribute to the 
observed phenomena as those formed before were expo- 
nentially diluted. Furthermore, whatever the considered 
formation mechanism, either through the usual collapse 
of high density-contrast primordial gaussian fluctuations 
or for near critical phenomena |4^], the PBH mass at 
the formation epoch is close to the horizon mass at 
the same time. It cannot be larger as the considered 
points would not be in causal contact and it cannot 
be much smaller as they would, in this case, more 
probably have formed before (as taken into account in 
the usual Press-Schechter formalism). It means that 
if the evaporation process was detected, the Hubble 
mass at the reheating time should be small enough not 
to induce a cutoff in the PBH mass spectrum which 
would make the light black holes abundance totally 
negligible. The best upper limit available on the density 
of PBHs around = 5 x 10^'' g, taking into account 
both the details of the source term evolution and the 
background from galaxies and quasars, is currently: 
^PBHiM^) < 3.3 X 10-9 m. 

Fortunately, some hope for future detection is still pos- 
sible thanks to antideuterons: those nuclei are expected 
to be very rarely formed by spallation processes below a 
few GeV for kinematical reasons. The threshold for an 
antideuteron production is E ~ 17 nip (total energy) in 
the laboratory, 2.4 times higher than for antiproton pro- 
duction. The center of mass is, therefore, moving fast 
and it is very unlikely to produce an antideuteron nearly 
at rest (in the 100 MeV - 1 GeV range) in the laboratory. 
On the other hand, they could be emitted in this energy 
range by evaporating PBHs and could be probed by the 
new generation of cosmic-ray detectors: the AMS exper- 
iment 01 and the GAPS projectE^. To obtain this 
result, a coalescence model (see for a review) was 
used, based mainly on phase space considerations: the 
antideuteron density in momentum space is proportional 
to the product of the proton density with the probability 
of flnding a neutron within a small sphere of radius po 
around the proton momcntmn. Thus: 



47r . / d^Np 



d^Nr, 



where po is the coalescence momentum whose uncertainty 
window is of the order of 60-280 Mev in extreme cases. 
The Hawking spectrum has then been convolved with the 
fragmentation functions, as obtained with the PYTHIA 
Monte-Carlo simulation of the Lund model: 
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where dgji^{Q, E,po)/dE is the number of antideuterons 
formed with an energy between E and E + dE by a 
partonic jet of type j and energy Q, evaluated with the 
coalescence model for a given momentum pq, aj is the 
number of degrees of freedom, s is the spin, and is 
the absorption probability. This coalescence condition 
(finding an antiproton and an antineutron within the 
same jet with a momentum difference smaller than po) 
was directly tested in the p — n center of mass frame 
as Po is not Lorentz-invariant and implemented within 
the PYTHIA simulation. This individual flux in then 
convolved with the PBH mass spectrum. To obtain the 
top of the atmosphere (experimentally measurable) spec- 
trum, the emitted antideuterons have been propagated 
within the Galaxy using the diffusion model of [s^ . 
briefly recalled in the appendix at the end of this paper. 
Finally, the resulting flux were solar-modulated in the 
force-field approximation. 

Figure |21 shows the possible values of the reheating 
temperature as a function of the density of PBHs at 
5 X 10^'* g for different PBH-induced antideuteron flux 
at 100 MeV (ranging from 2 x lO""^ m-^s-isr-^GeVS 
the maximum value consistent with the gamma-ray 
upper limit, down to 2 x 10~^° m~^s~^sr' 
They were obtained with conservative values of all 
the free parameters entering the model, astrophysical 
quantities being totally bounded by an exhaustive study 
of the heavy nuclei data As expected, there is a 

degeneracy between the D flux and flpBH- the same 
amount of particles can be produced either by a high 
normalization of the black hole spectrum and a cutoff in 
the high mass range {i.e. a low reheating temperature 
value) of by a low normalization of the black hole 
spectrum and a cutoff in the low mass range {i.e. a 
high reheating temperature value). This means that, 
in the case of detection, it should be possible to give a 
lower limit on the reheating temperature. Of course, the 
larger the antideuteron flux, the better the constraint 
on TjiH. As shown on this figure, for a fixed value of 
the fiux, whatever the value of flpBH, the reheating 
temperature value cannot be arbitrarily low since the 
mass spectrum cannot be cut much above masses 
roughly corresponding to temperatures of the order of 
the D mass {i.e. Tbh ^ a few GeV and T^h ~ a few 
10^ GeV). The other way round, whatever the value of 
TnH, the density of black holes cannot be arbitrarily 
low since even without any cutoff the source term must 
remain high enough to account for the considered flux. 
Naturally, this approach assumes that the measured 
antideuterons are indeed produced by evaporating black 
holes. The only other serious candidate as a source of 
light antinuclci in the low energy range are annihilating 
supersymmetric particles. It has been demonstrated |43 
that only neutralinos with masses around 100-200 GeV 
could contribute to the observed antideuteron flux. As 
this mass range will be probed by the Large Hadron 
Collider, it should be possible to distinguish between 
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FIG. 4: Lower limit on the gravitino mass as a function of 
the measured antideuteron flux for three different branching 
ratios. 



antideuterons induced by PBHs and by SUSY particles 
(some reconstruction problems could occur if the mass 
spectrum is strongly degenerated, especially between the 
lightest neutralinos and charginos, but this would hide 
the lightest supersymmetric particles only for masses in 
the TeV range). 

In the case where they are indeed coming from black 
holes. Fig. 13 gives the reheating temperature value as a 
function of the measured D flux. This result was ob- 
tained by varying values of the 100 MeV antideuteron 
spectrum combined with the upper limit coming from 
and {VlpBH < 3 x 10"^) for the correspond- 
ing reheating scale (evaluated by the previously given 
method). As expected, the limit becomes more stringent 
when the measured flux is higher and diverges when it 
goes to the maximum allowed value (otherwise it contra- 
dicts previously given limits). When compared with the 
upper bound coming from Big-Bang nucleosynthesis, this 
translates into a lower limit on the gravitino mass 
This can be derived by solving the Boltzmann equation 
for the gravitino number density (211] : 



dn 



3/2 



dt 



< E. 



3/2 > "^3/2 



where H is the Hubble parameter, Urad — C(3)r'^/7r^ 
is the number density of the scalar bosons in thermal 
bath, Vrei is the relative velocity of the scattering 
radiation, m^^2/{< >) is the averaged Lorcntz 

factor, T^/2 is the lifetime of the gravitino (computed 
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from the supergravity lagrangian ^u\) and T,tot is the 
total cross section (computed in the MSSM framework). 
Gravitinos are then assumed to decay mostly into 
photinos and photons, whose pair scattering off the 
background radiation, photon-photon scattering, pair 
creation on nuclei, compton scattering, inverse compton 
scattering of e+/e^ and induced leptonic cascades are 
taken into account. Requiring that the subsequent 
photo-dissociation of light elements does not modify the 
Big Bang Nucleosynthesis scenario beyond experimental 
constraints, the upper limit of the reheating temperature 
can be numeri cally computed as a function of the 
gravitino mass Figure ^ gives this bound as a 

function of the measured antideuteron flux at 100 MeV 
for three different branching ratios B of gravitinos into 
photons and photinos ranging from 0.1 (lowest curve) 
to 1 (upper curve). As the reheating temperature lower 
limit is extremely sensitive to the gravitino mass in 
the 100 GeV - 1 TeV range the curves are quite 
flat, except when the required value of Trh enters the 
diverging region. Although the accurate value of B is 
model dependent, it can safely be taken as lying in the 
0.1-1 range, as usually assumed in most studies. Once 
again, if a "thermal- like" inflation phase occurred, those 
limits do not stand anymore but could lead to important 
indications in favor of such a scenario if the gravitino 
was independently shown to be lighter than those values. 

It is important to notice that a great amount of work 
has also been recently devoted to the non-thermal pro- 
duction of gravitinos and moduli fields (dilaton and mod- 
ulus fields appearing in the framework of superstring the- 
ories which acquire mass through the nonperturbative ef- 
fects of the supersymmetry breaking) . Most papers claim 
that the upper limit on the reheating temperature must 
be drastically decreased (by up to 7 orders of magnitude 
[sH'l. Those results being still controversial, they were 
not taken into account in this work but they can only 
reinforce our conclusions and improve our limits. 



IV. PROSPECTS AND CONCLUSION 

It must be pointed out that such possible constraints 
on the gravitino mass can be translated into constraints 
on more fundamental parameters, making them very 
valuable in the search for the allowed parameter space 
in grand unified models. As an example, in models lead- 
ing naturally to mass scales in the 10^-10'^ GeV range 
through a specific dilaton vacuum configuration in su- 
pergravity, the gravitino mass can be related with the 
GUT parameters 
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With Mgut ~ 10^^ GeV and a gauge coupling a^f/T ^ 
1/26. The superpotential value in the dilaton direction 
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FIG. 5: Lower limit on the coupling constant A as a function 
of the reheating temperature. 



defines the magnitude of the coupling constant A of the 
self- interacting 24 multiplet. Figure shows how the 
lower value on A evolves as a function of the reheating 
temperature which could be probed by the previously 
given method, for three different branching ratios. Al- 
though not very constraining, this lower limit of the or- 
der 1.4 X 10~^ over the full tested range for S = 1 could 
be one of the first experimental constraints on A. 

The next generation of CMB experiments will face 
a new situation. Important efforts are devoted to the 
search for the polarization B mode 0| and the sen- 
sitivity should reach scales of inflation of order 10^^ — 
10^6 GeV. This value is slightly higher than the GUT 
scale if supersymmetry is ignored (i.e if gravitinos pro- 
duction is expected not to have occurred), and slightly 
lower than the GUT scale if supersymmetry is taken into 
account {i.e. in the case gravitinos are expected to be 
produced by scattering processes). Considering that the 
grand unified scale is the highest natural value for the 
reheating temperature, this means that, if a significant 
amount of entropy was not released after the moduli pro- 
duction, it should not be possible to detect those tensor 
modes in both scenarios. 

On the other hand, cosmic-ray experiments could be 
sensitive enough to investigate the allowed reheating 
temperatures if small black holes were formed at the 
end of inflation. In this case, important limits could be 
derived on the gravitino mass and on the related GUT 
parameters. 
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V. APPENDIX : ANTIDEUTERON FLUX 
COMPUTATION 

In this two-zone approach, the geometry of the Milky 
Way is a cyhndrical box whose radial extension is i? = 20 
kpc from the galactic center, with a matter (stars) disk 
whose thickness is 2h — 200 pc and a diffusion halo whose 
extent is the major source of uncertainty (taken into ac- 
count in the analysis). The five parameters used are Kq, 
5 (describing the diffusion coefficient K{E) — Kq(3R^), 
the halo half height L, the convective velocity Vc and 
the Alfven velocity Va- They are varied within a given 
range determined by an exhaustive and systematic study 
of cosmic ray nuclei data The same parameters as 
employed to study the antiproton flux |37| | are used again 
in this analysis. The antideuterons density produced by 
evaporating PBHs per energy bin ipQ obeys the following 
diffusion equation: 



written as 



V, 



d 



K 



d 



ipQ{r,z,E)- 



2h5{z)Vr>^I,{r,Q,E) = qP"^{r,z,E) 

where qP""^{r, z, E) corresponds to the source term. The 
total collision rate is given by = nncrj^f^VQ where 
ajjH is the total antideuteron cross-section with protons 
and the hydrogen density, assumed to be constant all over 
the disk, has been fixed to nn = 1 cm~'^. 

Performing Bessel transforms, all the quantities can be 
expanded over the orthogonal set of Bessel functions of 
zeroth order: 



1=1 

and the solution of the equation for antideuterons can be 



N. 



D,prira 



(0) = exp 



V^{L) 



2K J AjSinh(5ii/2) 



where 



= 2 exp (^(L - z')) sinh (f (i - z')) q^"^{z')dz' 



y2 /■2 ^ 1/2 



Vc + K Si colli 



In this model, energy changes (predominantly ioniza- 
tion losses, adiabatic losses and diffusive reacceleration) 
are taken into account via a second order differential 



equation for N 



D ,prim 



The spatial distribution /(r, z) 



of PBHs was assumed to follow 



Hi 



Rl + r-^ + z2 

where the core radius Rc has been fixed to 3.5 kpc and 
Rq—8 kpc. This profile corresponds to the isothermal 
case with a spherical symmetry, the uncertainties on Rc 
and the consequences of a possible flatness have been 
shown to be irrelevant in [37 
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